Local stability implies global stability for the 2-dimensional Ricker map

Consider the difference equation $x_{k+1}=x_k e^{\alpha-x_{n-d}}$ where $\alpha$ is a positive parameter and d is a non-negative integer. The case d = 0 was introduced by W.E. Ricker in 1954. For the delayed version d >= 1 of the equation S. Levin and R. May conjectured in 1976 that local stability of the nontrivial equilibrium implies its global stability. Based on rigorous, computer aided calculations and analytical tools, we prove the conjecture for d = 1.Comment: for associated C++ program, mathematica worksheet and output, see http://www.math.u-szeged.hu/~krisztin/ricke

Ab-initio molecular dynamics simulation of hydrogen diffusion in α\alpha-iron

First-principles atomistic molecular dynamics simulation in the micro-canonical and canonical ensembles has been used to study the diffusion of interstitial hydrogen in $\alpha$-iron. Hydrogen to Iron ratios between $\theta=1/16 and 1/2 have been considered by locating interstitial hydrogen atoms at random positions in a$2 \times 2 \times 2$ supercell. We find that the average optimum absorption site and the barrier for diffusion depend on the concentration of interestitials. Iron Debye temperature decreases monotonically for increasing concentration of interstitial hydrogen, proving that iron-iron interatomic potential is significantly weakened in the presence of a large number of diffusing hydrogen atoms

Hydrogen site occupancy and strength of forces in nano-sized metal hydrides

The dipole force components in nano-sized metal hydrides are quantitatively determined with curvature and x-ray diffraction measurements. Ab-initio density functional theory is used to calculate the dipole components and the symmetry of the strain field. The hydrogen occupancy in a 100 nm thick V film is shown to be tetrahedral with a slight asymmetry at low concentration and a transition to octahedral occupancy is shown to take place at around 0.07 [H/V] at 360 K. When the thickness of the V layer is reduced to 3 nm and biaxially strained, in a Fe_0.5V_0.5/V superlattice, the hydrogen unequivocally occupies octahedral z-like sites, even at and below concentrations of 0.02 [H/V]

Self-organized criticality induced by quenched disorder: experiments on flux avalanches in NbHx_x films

We present an experimental study of the influence of quenched disorder on the distribution of flux avalanches in type-II superconductors. In the presence of much quenched disorder, the avalanche sizes are power-law distributed and show finite size scaling, as expected from self-organized criticality (SOC). Furthermore, the shape of the avalanches is observed to be fractal. In the absence of quenched disorder, a preferred size of avalanches is observed and avalanches are smooth. These observations indicate that a certain minimum amount of disorder is necessary for SOC behavior. We relate these findings to the appearance or non-appearance of SOC in other experimental systems, particularly piles of sand.Comment: 4 pages, 4 figure

Zero-field and Larmor spinor precessions in a neutron polarimeter experiment

We present a neutron polarimetric experiment where two kinds of spinor precessions are observed: one is induced by different total energy of neutrons (zero-field precession) and the other is induced by a stationary guide field (Larmor precession). A characteristic of the former is the dependence of the energy-difference, which is in practice tuned by the frequency of the interacting oscillating magnetic field. In contrast the latter completely depends on the strength of the guide field, namely Larmor frequency. Our neutron-polarimetric experiment exhibits individual tuning as well as specific properties of each spinor precession, which assures the use of both spin precessions for multi-entangled spinor manipulation.Comment: 12 pages, 4 figure

Positive approximations of the inverse of fractional powers of SPD M-matrices

This study is motivated by the recent development in the fractional calculus and its applications. During last few years, several different techniques are proposed to localize the nonlocal fractional diffusion operator. They are based on transformation of the original problem to a local elliptic or pseudoparabolic problem, or to an integral representation of the solution, thus increasing the dimension of the computational domain. More recently, an alternative approach aimed at reducing the computational complexity was developed. The linear algebraic system $\cal A^\alpha \bf u=\bf f$, $0< \alpha <1$ is considered, where $\cal A$ is a properly normalized (scalded) symmetric and positive definite matrix obtained from finite element or finite difference approximation of second order elliptic problems in $\Omega\subset\mathbb{R}^d$, $d=1,2,3$. The method is based on best uniform rational approximations (BURA) of the function $t^{\beta-\alpha}$ for $0 < t \le 1$ and natural $\beta$. The maximum principles are among the major qualitative properties of linear elliptic operators/PDEs. In many studies and applications, it is important that such properties are preserved by the selected numerical solution method. In this paper we present and analyze the properties of positive approximations of $\cal A^{-\alpha}$ obtained by the BURA technique. Sufficient conditions for positiveness are proven, complemented by sharp error estimates. The theoretical results are supported by representative numerical tests

A note on the convergence of parametrised non-resonant invariant manifolds

Truncated Taylor series representations of invariant manifolds are abundant in numerical computations. We present an aposteriori method to compute the convergence radii and error estimates of analytic parametrisations of non-resonant local invariant manifolds of a saddle of an analytic vector field, from such a truncated series. This enables us to obtain local enclosures, as well as existence results, for the invariant manifolds

The High-Flux Backscattering Spectrometer at the NIST Center for Neutron Research

We describe the design and current performance of the high-flux backscattering spectrometer located at the NIST Center for Neutron Research. The design incorporates several state-of-the-art neutron optical devices to achieve the highest flux on sample possible while maintaining an energy resolution of less than 1mueV. Foremost among these is a novel phase-space transformation chopper that significantly reduces the mismatch between the beam divergences of the primary and secondary parts of the instrument. This resolves a long-standing problem of backscattering spectrometers, and produces a relative gain in neutron flux of 4.2. A high-speed Doppler-driven monochromator system has been built that is capable of achieving energy transfers of up to +-50mueV, thereby extending the dynamic range of this type of spectrometer by more than a factor of two over that of other reactor-based backscattering instruments

Computational Complexity of Iterated Maps on the Interval (Extended Abstract)

The exact computation of orbits of discrete dynamical systems on the interval is considered. Therefore, a multiple-precision floating point approach based on error analysis is chosen and a general algorithm is presented. The correctness of the algorithm is shown and the computational complexity is analyzed. As a main result, the computational complexity measure considered here is related to the Ljapunow exponent of the dynamical system under consideration

The Complexity of Flat Freeze LTL

We consider the model-checking problem for freeze LTL on one-counter automata (OCAs). Freeze LTL extends LTL with the freeze quantifier, which allows one to store different counter values of a run in registers so that they can be compared with one another. As the model-checking problem is undecidable in general, we focus on the flat fragment of freeze LTL, in which the usage of the freeze quantifier is restricted. Recently, Lechner et al. showed that model checking for flat freeze LTL on OCAs with binary encoding of counter updates is decidable and in 2NEXPTIME. In this paper, we prove that the problem is, in fact, NEXPTIME-complete no matter whether counter updates are encoded in unary or binary. Like Lechner et al., we rely on a reduction to the reachability problem in OCAs with parameterized tests (OCAPs). The new aspect is that we simulate OCAPs by alternating two-way automata over words. This implies an exponential upper bound on the parameter values that we exploit towards an NP algorithm for reachability in OCAPs with unary updates. We obtain our main result as a corollary